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MALMHEDEN'S THEOREM REVISITED 
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Abstract. In 1934 H. Malmhcdcn [15] discovered an elegant geo- 
metric algorithm for solving the Dirichlet problem in a ball. Al- 
though his result was rediscovered independently by Duffin [7] 23 
years later, it still does not seem to be widely known. In this pa- 
per we return to Malmheden's theorem, give an alternative proof 
of the result that allows generalization to polyharmonic functions 
and, also, discuss applications of his theorem to geometric proper- 
ties of harmonic measures in balls in M" . 

1. Introduction 

In 1934, H. Malmheden, a doctoral student of M. Riesz in Lund, 
proved that the value of a harmonic function u at any point P in a 
disk, or in a ball in M"-, n > 3, can be computed from its values on 
the boundary according to the following algorithm. Take an arbitrary 
chord L through P; calculate the value at P of the linear function i 
on L that interpolates the values of u at the endpoints of chord L; 
and, finally, calculate the average of the values i{P) over all chords 
L through P (cf. Theorem 12.11 for the precise statement), [15], also 
see [2-6]. 

To fix the ideas, if we denote, say in M^, by A{6) the value at P ob- 
tained by linearly interpolating the boundary values of a harmonic func- 
tion u in the disk at the endpoints of the chord through P making an 
angle 6 with the positive direction of the x-axis, Malmheden's theorem 
yields that u{P) = (l/27r) Jq"" A{e) d9. It follows then from Malmhe- 
den's theorem that u{P) < m.axA{9). In [3,4] J. Barta attempted to 

extend the latter inequality to all convex, or even star-shaped regions. 
In [17], Weinberger showed that no such inequality, even allowing a 
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multiplicative factor on the right, can hold at all points of a convex 
region unless the region is a disk. The analogous assertion is true in 
higher dimensions. Thus, the converse to Malmheden's theorem also 
holds; i.e., the assertion of his theorem is correct only if the domain is 
a ball. (Note, in passing, that the complex-analytic version of Malmhe- 
den's theorem is well-known under the name of Bochner - Martinelli 
formula. It provides holomorphic extensions of analytic functions from 
the boundaries of arbitrary smoothly bounded domains in and is 
obtained by averaging one- dimensional Cauchy integral representations 
-cf., e.g., [10].) 

With this paper we hope to return the beautiful result of Malmheden 
to mathematical folklore since it appears that the result is not widely 
known. Moreover, we present an alternative approach to Malmheden's 
theorem that allows to generalize the result to the polyharmonic func- 
tions, cf. [2, 5, 6]. We also present applications of Malmheden's theorem 
yielding nice geometric properties of the harmonic measure in balls in 
M". 

Regarding the history of the result, it must be mentioned that Duf- 
fin [7] has re-discovered Malmheden's theorem more than two decades 
later. His proof, although obtained independently, is essentially that 
of [15]. 

The paper is organized as follows. In ^and ^we give two different 
proofs of the Malmheden theorem and its generalization to fc— dimensional 
cross-sections of the ball passing through a given point P. The first 
proof in ^is essentially Malmheden's original proof. The second proof 
seems to be new and is the crux for further generalizations of Malmhe- 
den's procedure to polyharmonic functions. In §l]we show the converse, 
i. e., that balls are the only domains in for which Malmheden's the- 
orem holds. In [|5]we discuss Malmheden's theorem in relation to some 
geometric properties of harmonic measure in the ball. In ^we extend 
Malmheden's theorem to polyharmonic functions in M". In 2 dimen- 
sions this was done by Barta, cf. [2,5,6]; his proof is quite different 
from ours. A different converse theorem, assuming that Malmheden's 
algorithm only reproduces harmonic functions at one point and, under 
additional condition of central symmetry, implying that the domain is 
a ball is given in §71 We finish with some additional remarks in §Hi 

2. A PROOF OF Malmheden's theorem 

Let us begin in a more general setting. Let Q he a. convex open set 
in R" with boundary F. Let / be a continuous function on F, P G Q. 
Draw a chord L through P intersecting F in exactly two points Qi, Q2- 
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Denote by /i, /2 the values of / at Qi, Q2, respectively. Interpolate 
the values /i, /2 by the unique linear function on L that we denote by 
i. Clearly, the value i{P) is a continuous function of P and L. Holding 
P fixed let u{P) denote the average of i{P) over all lines L. u{P) is 
a continuous function in Q. Note, that if / is a restriction of a linear 
function onto F, then f = u in Q. 

Clearly, u{P) ^ / (Qo) when P Qo, Qo e F. 

Theorem 2.1 (Malmheden — [7,15]). If n is a ball in M", then u is 
harmonic in Vl, and hence solves the Dirichlet problem (for the Lapla- 
cian) in Q with data f . 

Proof. In view of the preceding remarks it remains to show that u is 
harmonic. Let L, P, Qi, Q2 be as above, ri = |PQi|, r2 = |-PQ2|- 
Then 

(2.1) ^{P) = !ili±lili, 

ri + r2 

First, assume n = 2. Let 6 be the (polar) angle that the chord L 
makes with the positive direction of the x-axis, so ri{6 + tt) = r2{6) 
and fi{9 + vr) = f2{0)- Then, we have, 

2tt 



uiP) ^ I J HP) 

(2.2) 



de 





1 f r^f: 



TT J ri + r2 



— d9. 







Let if be the angle between the (inner) normal n to F at Qi and the 
chord L. Then, comparing infinitesimal arclengths on F and the circle 
of radius ri, centered at P, we observe that rid9 = cos ipds, where ds 
is the arclength on F. Thus, from (12.21) it follows that 

(2.3) u(P)='-h^f,ds-l[^nds, 

n J ri n J ri+r2 



where c is the length of F. 

cos (y9 

The first term in (12.31) is harmonic since is nothing else but 

n 

d 

log|P — Qi\, so the first integral is the "double layer" potential 



of /i (cf. [12]). Finally, we use the hypothesis that F is a circle of 
radius R. Then vi +r2 = 2Rcosip, and the second integral is simply a 
constant, hence u is harmonic in Q. 
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In higher dimensions, the argument is almost identicaL For the sake 
of clarity of notation, we take n = 3. Then, fl2.2p and (I2.3p become 

2tt J ri + ra 



r 



and, hence. 



Z-K J r{ 2tc J ri + r2 ri 

r r 

Here, dS is the surface area measure on F and d6 is the solid angle 
measure at P. Once again, the first integral is the double layer potential 
of /i, while when is a ball of radius R, ri + r2 = 2Rcosip, so the 
second integral becomes a "single-layer" potential of /i, harmonic in 

n. 

(Of course, (12 .Sip and (12.31) correspond to very well known repre- 
sentations of solutions of the Dirichlet problem in the ball as sums of 
potentials of double and single layers in dimensions 3 and higher, and 
as (up to an additive constant) a potential of a double layer in the disk 
(cf. [12] or [13]). □ 

Let us sketch here a different proof of Theorem 12.11 that will allow 
us in ^to generalize Malmheden's algorithm to biharmonic functions. 

2nd Proof of Theorem 12. 1[ It is more convenient to translate our coor- 
dinate system so the origin is now at P while f2 = {xGM"': |a; — c| <1}, 
where c : |c| < 1 is the center of the ball. We wish to show that for any 
homogeneous harmonic polynomial H in dQ, Malmheden's procedure 
produces the number H{0). Clearly, this will imply the result since 
the restrictions of harmonic polynomials to the sphere dQ are dense in 
C{dfl), while linearity of Malmheden's algorithm yields the result for 
all harmonic polynomials as long as it holds for homogeneous polyno- 
mials. Finally, since Malmheden's theorem obviously holds for linear 
polynomials we may assume that m := degH > 1. and need to verify 
that the algorithm produces the number = H{P). Fix e a unit vector 
in M", and let L = L^he the line through the origin parallel to e. Let 
t be the running coordinate on L^. meets dQ at two points Qi, Q2 
where t is equal to the roots of the equation 

(2.4) |te-c|^ = l, _ 2 < c,e > t + (|c|^ - 1) = 0. 

Along Le, H (te) = t'^H{e). If we denote by a, 6 G M the values of 
t corresponding to Qi, Q2 and determined by (12. 4p . then the linear 
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interpolant of H along the segment of Le between Qi, Q2 equals 



and is equal at t = to 



Now, we note from (12 .4^ that a6 = |cp — 1, independent of the chosen 
vector e. (Of course, this is a well-known theorem from the euclidean 
geometry). So the cofactor of H (e) in (12.61) is a symmetric polynomial 
of degree m — 2 of the roots a, b of (12. 4p . Hence, by a standard result 
on symmetric functions, it is a polynomial of degree at most m — 2 
of the coefficients of the quadratic (12.41) . Now observe that the only 
coefficient dependent on e in (12.41) is that of t and it is of degree 1 in 
the coordinates ej of e. Summarizing, the cofactor of H (e) in (12.61) is 
a polynomial in ej of degree at most m — 2. Hence, the cofactor as a 
function of e is orthogonal to the homogeneous harmonic polynomial 
H (e) over the unit sphere, thus Malmheden's procedure for H indeed 
produces = H{0). The last assertion follows from the fact that a 
polynomial of any degree k (in our case, k = m — 2) matches on the 
sphere a harmonic polynomial of degree < k, which, of course, is the 
sum of homogeneous harmonic polynomials. Homogeneous harmonic 
polynomials of different degrees are orthogonal on the unit sphere, cf. 
e.g., [9,12,14]. The second proof is now complete. □ 

3. An EXTENSION OF Malmheden's theorem to 

CROSS-SECTIONS OF HIGHER DIMENSION 

Theorem 12.11 has an almost immediate generalization if one replaces 
chords through a point P by fc-dimensional cross-sections. (For the 
notational convenience we shall often represent P in the coordinate 
form as a; = (xi, . . . , x„)). More precisely, let Q he a convex bounded 
domain in R" with smooth boundary F. Let P{x,y), x & Q, y & 
F, denote the Poisson kernel for Q, so any harmonic function in Q, 
continuous in Q, can be represented by its boundary values f{y) on F 
as 



(2.5) 




(3.1) 




Let, as usual, the Green function in f2 be 



(3.2) g{x,y) = \ if 
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where Ux{y) is harmonic in Q, and g{x, y) \y&T= 0. oJn denotes the area 
of the unit sphere in M". Then, as is well-known, 

d 

(3.3) P{x,y) = —g{x,y), 

when Uy is the inner normal to F at ?/ G F. For Q = {\x\ < 1}, 



(3.4) Pix,y) 



1 1 - |x 



2 



\x - y\ 



(cf. [12], e.g., for more details). 

For n = 1 there is a formula playing the role of (13. ip and (13. 4p , where 
the "domain" is the interval [—1, 1]. Then 



x) [x + 1, 



i[/(l)(l + x) + /(-l)(l-x)] 



is the linear function matching given data / at ±1. 

Fix k: 1 < k < n and consider a fc-dimensional plane a through P. 
f^Q, := n n a is a convex domain in M^. Let P" denote the Poisson 
kernel for Q^- Replacing in (13.11) P{x, y) by P°, f hj fa = f \anr, and 
dS{y) by dSa, the Lebesgue measure on aflF, yields the solution of the 
Dirichlet problem P"{fa) for Qa with the data fa with respect to the 
fc-dimensional Laplacian. Fix x E fl and consider the Grassmann man- 
ifold G{k, n) of all fc-dimensional planes a through x. The orthogonal 
group SO{n) acts naturally on G{k, n) and the invariant Haar measure 
on SO{n) induces the unique normalized measure dmk{a) on G{k,a). 
mfc(a) is invariant with respect to all rotations of R"' with the center 

X. 

Now keeping the notations from §21 and in view of the note following 
(13.41) . we can restate Theorem 12.11 and rewrite (13. ip . Ifi7 = {|x|<l}is 
the unit ball, then any harmonic function u in Q, continuous in Q with 
u |r= / can be represented via solutions of one- dimensional Dirichlet 
problems on 1-dimensional lines L passing through x & fl. Namely, 
denoting by fi := / Unr, we have: 

(3.5) u{x)= J P''{fL)dm,{L). 

LeG^(l,n) 

This yields an immediate corollary. 

Theorem 3.1. Let n = {\x\<l} he the unit hall, f e C(T), T = dn is 
the unit sphere. The harmonic extension u of f to Q can he represented 
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in terms of harmonic extensions of f into k- dimensional sections ofVt: 
(3.6) u{x) = j P° (falrna) ix)dmk{a). 

Proof. From (13.51) and Fubini's theorem we get 



u[x) = 

LGG^(l,n) 



j P^(/l IrnL)f^mi(L) 
j J if L Irni) dmi{L) dnikia) 



= J P''ifa\Tna)dmj,{a). □ 



a€Gx{k,n) 



□ 



4. The converse 



Theorem 4.1. If Q is a convex domain and the Malmheden procedure 
described in ^ produces for each continuous function f the solution 
of the Dirichlet problem with data f , then Q must be a ball. 

For the sake of clarity, we restrict ourselves to the cases n = 2, 3. 
The extension to higher dimensions is straightforward. 

We shall use the notation from ^ The hypothesis implies that the 
function u{P) given by ([23D or (HI is a harmonic function of P. Note 



that the first term, as we've observed in §2], is ALWAYS harmonic in 
Q since it is equal to the double layer potential. Thus the hypothesis 
yields that the second term is harmonic as a function of point P E Q 
as well. Therefore the functions 

4.1 1 — 11 ^ 

ri + r2 ri(ri+r2) 

are harmonic (as functions of P) in Q. 

Put the origin at the point Qi. Let P = rit, where t : |t| = 1 is a 
point on the unit sphere centered at Qi. 

Let us represent F := dfl near Q2 by its "polar" equation: p = p(t), 
P = \QiQ2\- Then, (14. ip yields that functions 

(4.2) (i) (ii) i 

p{t) ri pit) 

where (, ) denotes the scalar product and uq^ is the (inner) normal to 
F at Qi, are harmonic functions of ri and t. 
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In (i), the function depends on the polar angle t only, and being 
harmonic, forces it to be linear. Yet, it must be single- valued, thus is 

a constant. Denote it by — . Thus, p{t) = 2i?cos (((t, ngj)) which is 

an equation of the circle. 

For (ii), a simple calculation, or a quick check with [12, p. 141, Ex. 4] 
const const 

yields that ( — n > 3, m general) are the only homogeneous 

ri r'l 

harmonic functions of degree —1 (2 — n, respectively). Hence, again, 

from (14. 2p . (ii) we infer that ^ ' .'^^'^ = const, i.e., F is a sphere. □ 

pit) 

Remark 4.2. Weinberger [17] actually showed a stronger converse. To 
fix the ideas, we shall describe his result for n = 2. 

We keep the same notation as above. If Malmheden's procedure 
produces the solution of the Dirichlet problem, then putting the origin 
at point P E Q, assuming Q to be convex, and the boundary F to be 
given by F := {r = r{6), < 6 < 27r}, we can rewrite (12. 2p from (12.11) 
as 

U 3) u(P) = - (rje + n),e + n) + rie + 7r)^(r(g), 6) 



Equation (14.31) implies 

r{9)u{r{9 + vr), ^ + vr) + r(0 + 'n)u{r{e), 6) 



(4.4) u{P) < max ■ 



r{9) + r{9 



Weinberger showed that if a weaker version of (14.41) with a constant 
factor on the right hand side holds for any P E Q, Q being a convex 
domain, then f2 is a disk. This result refuted the conjecture of Barta 
[4, 5] that (14.41) should hold in general convex domains. 

Remark 4.3. Going over the proof of Theorem 14. H one immediately 
notes that the hypothesis can be somewhat weakened: 

(i) Instead of convexity, it suffices to assume that the domain Q is 
star-shaped with respect to all points Q in a neighborhood of 
the fixed boundary point Qi G dQ. Equivalently, that placing 
the origin at any such Q, the (smooth, say C^) boundary dQ is 
given by the "polar" equation p = p(t), where t : \t\ = 1 runs 
over the upper hemisphere of the unit sphere. 

(ii) Instead of assuming that Malmheden's algorithm successfully 
solves the Dirichlet problem for any data /, it suffices to assume 
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that merely for any continuous data supported in a neighbor- 
hood of Qi. We omit the details. 

5. Discussing Malmheden's Theorem and Harmonic 

Measure 

We keep the same notation as in [|2l If f2 is a convex domain in i?" 
and P is an interior point, the metric ratio Rp associated with the pair 
P) is the function defined on d^l by: Rp{Qi) = j^^, where Q2 
denote the second point of intersection with dQ of the line L joining P 
to Qi. As in §21 we note that Rp is a continuous function on dQ with 
values strictly between and 1. It is constant (i.e., identically equal 
to 1/2) if, and only if, P is a center of symmetry of Q (i.e., all chords 
through P are bisected there). Let wp denote the harmonic measure 
on dfl evaluated at P and Ap denotes the subtended angle measure on 
dQ with respect to P. Malmheden's theorem then can be reformulated 
as follows. 

Theorem 5.1. If Q is a ball, then for any P G f2 , Rp is the Radon - 
Nikodym derivative of the harmonic measure wp with respect to Ap. 

For the proof one only needs apply Theorem 12. II to the characteristic 
functions of arbitrary (say, relatively open) subsets of dfl. 

Thus, Malmheden's theorem can be seen as a procedure to (when Q is 
a ball ) compute harmonic measure from "purely geometric" quantities, 
lengths and (solid) angles, by integration. Note that the corresponding 
derivative when is replaced by surface measure on dfl is given by 
C{P)/\PQi\'^, the Poisson kernel ([33]) evaluated at P and Qi. For 
the unit ball centered at the origin O, C{P) = ^(1 — lOP]"^), where 
Un is the surface area of the unit sphere - cf. [9, 12] . Modulo the 
known relation between measures Ay^ and the surface area measure 
Malmheden's theorem is thus equivalent to Poisson's formula for the 
solution to Dirichlet's problem for the ball. 

Against this background it is interesting to note the following fur- 
ther connections between harmonic and subtended angle measures in 
the ball, mediated by Malmheden's theorem. It implies other elegant 
properties of harmonic measure which, in dimensions > 2, are rarely 
noted. 

Corollary 5.1. Let P be a point in the unit ball fl := {x : |x| < 1}, 

and let a double cone Kp with vertex at P cut out "spherical caps" U 
and V from the unit sphere S. Let Wp{E) denote the harmonic measure 
of a set E (Z S evaluated at P. Then, wp{U) + wpiV) equals twice the 
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(normalized) solid angle at the vertex P of K. (The normalizing factor 
equals —, where a;„ denotes the surface area of the unit sphere in MJ^.) 



Keeping the same notation as in Corollary 15.11 we also have the 
following. 

Corollary 5.2. // we consider a system of masses consisting of caps 
U, V, each endowed with the harmonic measure wp{U), wp{V), re- 
spectively, then the center of mass is at P. 

Remark 5.2. Corollary 15. II is well known in 2 dimensions — cf. [16, Ch. 
IV, §2] and can be used to give an even shorter proof of Malmheden's 
theorem. Corollary 15. 2[ even in the two dimensional case, seems not to 
have been noticed before. 

Proof of Corollaru \5.1[ Fix P ^ Q. Let J = Jp be the involution of S 
that maps a point x of the unit sphere S onto the second point where 
the line Px meets S. Let / G C{S) be self-involutory with respect to 
J, i.e. f{x) = f{Jx), let us denote by dAp the normalized sohd angle 
measure at the vertex P. Then, Theorem 12.11 implies that 

(5.1) f{x)dwp{x) = I f{x)dAp{x). 



Thus, (15. ip holds for all self-involutory / G C{S). But the character- 
istic function Xuuv of the union UUV is obviously a bounded pointwise 
limit of the self-involutory (w.r.t. J) functions in C{S). Hence, by the 
bounded convergence theorem (15. ip also holds for Xuuv- This is the 
conclusion of Corollary 15. 1[ □ 

Proof of Corollary \5.2[ If i{x) is a linear polynomial and / G C{S) : 
f{x) = f{Jx) is self-involutory with respect to Jp, then the Malmheden 
algorithm applied to {£f){x) (of course, w.r.t. P) produces 

(5.2) e{P) [ f{x)dApix). 



(This is seen at once since Malmheden's algorithm preserves linear 
functions while f{x) = f (Jpx)). 
By Theorem 12.11 we have then 

(5.3) £{P) J f{x)dAp{x)=i{P) J f{x)dwp{x). 

s s 
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So, from (15. 2p . we therefore infer that 



i{x) ■ f{x) dwp{x) = i{P) J f{x) dwp{x). 
s s 
Hence, for all linear functions i, we have 

/ ■ fix)) dwp{x) 

(5.4) i{P) = ' , . 

j f{x)dwp{x) 

s 

Substituting for i the coordinate functions Xj,j = l,...,n, we ob- 
tain from (15 .4^ at once that the center of mass of the mass density 
f{x) dwp{x) is at P. Taking for / the characteristic function Xuuv of 
the union U UV proves the corollary. □ 

Remark 5.3. (i) Corollary 15. ll has an independent proof in 2 dimen- 
sions. If P, with a slight abuse of notation, denotes a complex 
number in D, D = {|2;| < 1}, it is straightforward to calculate 

P - z 

the involution Jp{z) to be Jp{z) = Consider an arc 

1 — Pz 

U on the unit circle T. Its harmonic measure at P equals the 
harmonic measure of V = Jp{U) evaluated at (by conformal 
invariance of the harmonic measure, since Jp{P) = 0). The 
harmonic measure of an arc V evaluated at the origin equals 
the central angle associated with V normalized by the factor 
1/2tt. Thus, the harmonic measure wp{U) at P equals 

(5.5) wp{U) = ^, 

Zn 

where \V\ denotes the length of the arc V = Jp{U). Hence, 
recalling that JpiV) = U, we conclude that 

wpiU) + wpiV) = ^i\U\ + \V\) 

(5.6) = — [Ap := angle subtended by U w.r.t. P) 

TT 

= — {Ap := angle subtended by V w.r.t. P) . 

TC 

(The last two equalities are, of course, corollaries of Euclid's 
theorem about angles formed by the two chords through a point 
in the disk which states that their radian measure equals to 
the average of the radian measures of the subtending arcs of 
the circle.) This gives a different proof of Corollary 15.11 in 2 
dimensions. 
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(ii) From (i), one can easily deduce Theorem 12.11 for n = 2. In- 
deed, Corollary 15.11 yields that for every w G D with, as before, 
denoting the (normalized) subtended angle measure as seen 
from w, and for every analytic polynomial P{z)^ we have 

(5.7) j P{OdA^{O = liP{0) + P{w)). 

T 

w — z 

(Recall that Jw{z) = —, J«,(0) = w, Jwiw) = 0.) Ob- 

1 — wz 

viously, (15. 7p then holds for all functions / analytic in D and 
continuous in ©. Let / be any such function. Parameterize 
the lines through w by their intersection points z with the unit 
circle T. It is easy then by solving a system of linear equations 
to find for the line i := i{z) the linear interpolant for the values 
f{z), f {Jw{z)). At w it equals 

(5.8) nj /(^)-/('^-(^)) + ^fiJM)-fi^)JM _ 

z-J^{z) z-J^{z) 

Clearly for a fixed w G D and / analytic, say, in a neighbor- 
hood of D (15.81) produces an analytic function F{z). Remem- 
bering that maps to w and w to 0, one easily calculates 
that -F(O) = F{w) = f{w). Hence, applying Malmheden's pro- 
cedure to /(2;), z G T and then using F{z) from ( 15.80 and (15.70 
we obtain the value f{w). Thus, Malmheden's algorithm ap- 
plied to, say, any polynomial in z, produces the "correct" value 
at w. 

Separating real and imaginary part yields Malmheden's the- 
orem in the disk, i.e., for n = 2. 

The following proposition can be viewed as a converse of Corollaries 
15.11 and 15.21 combined. 

Proposition 5.1. Let Q be a bounded convex domain in M" and P G fi. 

Assume that the harmonic measure wp = w{-,dQ, P) on T := dQ at 
P satisfies the conclusions of Corollaries 15.11 and 15. 2[ Then, P is a 
"Malmheden" point for Q , i.e., Malmheden's algorithm (cf. §§i and 2) 
applied to any continuous function f G C(T) on T produces the value 
at P of the solution to the Dirichlet problem in Q with data f onV . 

Proof. Let /(F) C C(F) consist of all self-involutory functions, i.e., such 
that / G C(F) : f{J{x)) = f{x), x G F. (Here, J = Jp is as above, the 
involution of F induced by P.) Denote by i? C C(T) the subset of / G 
C(F) for which the Malmheden algorithm produces the desired value at 



MALMHEDEN'S THEOREM REVISITED 



13 



P of the solution of the Dirichlet problem with data /. Our hypothesis 
implies that R D {/(r)}U{(£ ■/),/€ I{T),i is a linear function of Xi, . 
We need to show that R contains all polynomials, i.e., is dense in C(T). 
Since R is obviously closed the proposition will then follow. 

Fix a linear function i{x). Denote = i o J. We have i"^ = £ ■ 
+ — i ■ and, since £ + £ ■ G /(F), our hypothesis implies 
that for any g E I{n), f ■ g = i ■{£ + £*) g ~ {i ■£*) g E {/(F)} U 
{{e- f)J E 1} C R. So, {f-I{T)} C R. An induction argument 
shows that {i"^ ■ /(F)} C R for any integer m. Since the linear span 
of the set of powers of linear polynomials contains all polynomials the 
proposition follows. □ 

Remark 5.4. It seems a worthy question whether the two geometric 
properties of harmonic measure enunciated in Corollaries 15.11 and 15.21 
(i.e., the two hypotheses in Proposition 15.11) are actually independent 
of each other. We suspect they are but haven't been able to prove it. 



6. Malmheden's theorem for polyharmonic functions 

Let Q be the unit ball in M". F := dQ is the unit sphere, and let / : = 
f{x) be a smooth function, say even real- analytic, in a neighborhood 
of F. Then as is well-known (cf. [9]) the Dirichlet problem for the 
bi-harmonic operator 



(6.1) 



A^u = 0; 

u = f, Vu = V/ on F 



has a unique solution in Vt. Not going into technicalities, the reader 
may argue as follows. 

As is well-known — cf. [1], any function u : A^m = in a domain f2 
admits so-called Almansi expansion: 

(6.2) u = hi + \x\^h2, 

with hi and /12 harmonic functions in f2, uniquely defined by u. One 
can trivially adjust (16.21) when fi = {|x|<l}isa ball and rewrite (16.21) 

as 

(6.3) u = Hi + - 1) //2, 

//i,//2 being harmonic in Vt. Then, to solve the BVP (16.11) . we need 
to solve two consecutive Dirichlet problems: first, for Hi (= / on 
F), and then for H2: ^(/ — Hi) = 2H2, where ^ stands for the 
radial derivative. Uniqueness follows by application of Green's formula 
making use of the bi-harmonic Green's function — cf., e.g., [11]. 
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Now to formulate the analogue of Malmlieden's theorem for bihar- 
monic functions (i. e., those for which = 0) we proceed as follows. 

The one-dimensional analog of the bi-Laplacian is the operator (^^j , 
whose kernel consists of cubic polynomials. We can extend Malmhe- 
den's procedure for a fixed point P ^ Q to solutions of (16.11) . Again 
draw a chord L through P intersecting T = dQ at points Qi,Q2- Let 
t denote the real parameter along L, a and b being the values of t at 
Qi,Q2- Let C{t) denote the (unique) cubic polynomial such that the 
functional C(a), C'{a), C{b), C'{h) interpolate the corresponding val- 
ues of the data fhif) '■= f \l- (For derivatives of /l at a and b we 
take the values of the directional derivatives of the data / along L at 
those points.) Cl{P) is then the value of C at P. Holding P E Q 
fixed, let u{P) denote the average of Cl{P) over all lines L. u{P) is, of 
course, a continuous function in B. Also, if f{x) is a cubic polynomial, 
then u = f inn. Clearly, u{P) f (Qo), Vm(P) ^ V/ (Qo) when 

P ^ Qo, Qo e r. 

Theorem 6.1 ( [2,6]). Iffl is a ball in M", u is biharmonic in fl, and 
hence solves (16.11) . 

Proof. It is convenient to translate the coordinate system so f2 is a 
unit ball centered at some point c, while P is now the origin. By a 
standard approximation argument we may assume that the data f{x) 
is a polynomial since the latter are dense in the space C^(r) of con- 
tinuously differentiable functions on F. Moreover, we may also assume 
that f{x) is a biharmonic polynomial since the latter are also dense in 
the space of smooth biharmonic functions in the C^-metric. Lastly, by 
linearity of the Malmheden operator, we may assume that f{x) is in 
fact a homogeneous biharmonic polynomial of degree m > 3. 

Let us separate the following one-dimensional interpolation lemma. 

Lemma 6.1. Let a<0<b, m>3 and C{t) = Cm{t) be the (unique) 
cubic polynomial interpolating the values of t"^ and its first derivative 
at a and b. Then, Cm(0) is a homogeneous symmetric polynomial of 
degree m. Moreover, 

(6.4) CU0) = iabfqm-4ia,b), 

where Qm-i is a homogeneous polynomial of a and b of degree m — 4. 

Assuming the lemma and in view of our chain of reductions we need 
only show that for f{x) =: H{x), H{x) being a homogeneous bihar- 
monic polynomial of degree m > 4, the Malmheden algorithm produces 
the number = H{P), (Point P, recall, is at the origin). Let L, a, 
b be as described in the paragraph preceding Theorem 16. 1[ Note that 
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ab = constant that depends on P only as already noted in ^ (cf. the 
argument below (12.61) ) . From this, (16 ■4p and since if is a homogeneous 
polynomial it follows, by repeating the argument following (12.61) . that 
for each line L through P = defined by the directional vector e G F, 
Malmheden's procedure produces 

(6.5) const gm-4 (e) H (e) 

where the constant only depends on point P. 

Averaging (16. 5p over all chords C, i.e. over all directional unit vectors 
e G F produces as required. To verify this last assertion we simply 
note, since H is homogeneous, (16.21) implies that 

(6.6) H{x) = hi{x) + |xp/i2(x), 

where hi,h2 are homogeneous harmonic polynomials of degrees m and 
m — 2 respectively, qm-4, (as is true for any polynomial — cf. e.g., [14]), 
can be matched on the unit sphere F by a harmonic polynomial hq, 
deg hg < m—4: < m—2. From the well-known orthogonality of spherical 
harmonics (see, e.g., [14], [1], [9], [12]), we obtain our last assertion and, 
hence. Theorem 16.11 follows modulo Lemma [6.11 □ 

Proof of LemmaEH Let C^(t) := C{t) = At^ + Bt^ + Ct + D match 
the function together with its first derivatives at a and b, a < < b. 

Writing down four equations corresponding to these 4 interpolating 
conditions we see at once (by Cramer's Rule) that all the coefficients 
are rational functions of a, b of degree at most m + 6 symmetric with 
respect to a and b (interchanging a and b merely permutes equations 
in the system). 

Furthermore, solving this linear system of equations via Cramer's 
Rule we also observe that each coefficient A, . . . , D is a rational func- 
tion of a and b for which the degree of the numerator is at most m + 6, 
while the degree of the denominator which is the determinant of the 
system is precisely 6. Now these rational functions actually cannot 
have any finite poles for some complex values of a and b since Hermite 
interpolation polynomial Cm{t) exists and is unique for all complex 
values a, b. Thus, all the zeros of denominators in the rational expres- 
sions for the coefficients A, . . . , D of Cm{t) must cancel out. Hence, all 
the coefficients A, . . . ,D of Cm(t) are actually polynomials of degree 
(m-|-6) — 6 = m. Finally, since when a or b = 0, D = D{a, b) must have 
a double zero at the origin, it follows that D{a,b) = {ab)'^qm-4:{a,b), 
where qm-A is a symmetric polynomial in a, b of degree < m — 4. The 
lemma is proved and the proof of Theorem 16.11 is now complete. □ 
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Remark 6.2. 

(i) Obviously, the idea of the above argument originated from our 
second proof of Malmheden's original result — Theorem 12. II — 
given in §21 

(ii) Lemma 16.11 has a natural extension to higher order differential 



operators ( — ) , k > 3. Accordingly, with rather obvious 



modifications Theorem 16.11 extends to polyharmonic operators 



(iii) Theorem 13.11 also readily extends to polyharmonic operators. 

(iv) We do not know whether the converse to Theorem 16.11 similar 
to Theorem 14.11 also holds for polyharmonic functions. 

7. Another converse to Malmheden's theorem 

Let, as before, Q he a convex, bounded domain in M", F = dQ, 
P G fi. If Malmheden's procedure with respect to P as described in 
§2 applied to any, say, harmonic polynomial h yields the value h{P), 
it does not seem to imply that f2 is a ball. The problem, of course, is 
that it does not allow us to locate the center of the ball. However if 
one assumes that not only Malmheden's procedure applied to harmonic 
functions u in Q yields the correct value u{P), but also Q is centrally 
symmetric with respect to P, one can conclude that f2 is a ball centered 



The driving force for this result is the following extremely simple but 
useful observation regarding harmonic measures. 

Lemma 7.1. Let Qi ^ Q2 be two smoothly bounded, star-shaped do- 
mains in and assume that fii fl ^2 7^ and dfli fl dQ2 7^ 0- Let 
O & fli (1 fl2 be a point in f2i fl ^2- Let wi,W2 be harmonic measures 
on dQi, dQ2, respectively, evaluated at O. Let Ei = dQi \ Q2 be the 
portion of dQi that lies outside Q2, while F2 = dQ2 H ^1 is the por- 
tion of dQ2 that lies inside Qi. Similarly, define E2 = 8^2 \ ^i? and 

Fl = dninQ2- 




cf. [2,6]. 



at P. 



Then. 



(7.1) wi{Ei) <W2{F2), 
and, similarly, 

(7.2) W2 {E2) < wi (Fl) . 



(Both inequalities are strict.) 



Proof. We just prove (17.11) : the proof of (17. 2p is identical. 
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Let U = nin Q2, dU = FiU F2. On dU \ (Fi n F2) the (harmonic) 
functions Wi{Ei]x), W2 {F2]x) satisfy (by the maximum principle) the 
following: 

^^g^ Wi(Ei;x) |fi= 0; Wi{Ei]x) \f2<1 

W2 (F2;a;) |fi> 0; W2 (-F2; a;) |f2= 1- 

Thus, 

(7.4) W2{F2]x)>wi{Ei-x) 

on dU\ {Fi n F2). Since Fi r\F2 has measure zero on dU , and hence its 
harmonic measure is zero as well, the generalized maximum principle 
— cf., e.g., [9,12,16] — for bounded harmonic functions yields that 
(17.41) holds everywhere in U . This proves (17.11) . □ 

Theorem 7.1. Let VL he star-shaped domain in and O E Q be a 

point in fl. Assume also that Q is centrally symmetric w.r.t. O. If 
the Malmheden algorithm applied to O reproduces the values at O of all 
functions harmonic in Q and continuous in Q, then Q is a ball centered 
at O. 

Proof. First, observe that Lebesgue dominated convergence theorem 
immediately extends the hypothesis to all bounded harmonic functions 
whose boundary values are pointwise limits of continuous functions on 
dQ. Thus, in particular, the conclusion applies to harmonic measures of 
smoothly bounded open subsets of dfl (or, to subarcs in 2 dimensions). 
Let B be the ball centered at O with same volume as Q. If B = Q 
there is nothing to prove. Then, B ^ Q and since Vol(i?) = Vol(f2), 
dB n dQ 7^ 0. Since the Malmheden algorithm applies to Q at O, and 
Q is centrally symmetric with respect to O, the proof of Corollary 15.11 
implies that the harmonic measure on dfl evaluated at O is identical 
with the normalized solid angle measure subtended from O. Thus, 
applying Lemma 17.11 to the configuration Qi = Q, Q2 = B and O G 
Q n B, we arrive at the contradiction, since both harmonic measures 
on dfl and dB at O equal to the solid angle measure subtended from O 
and hence, must be the same for the sets Ei = dQ\B and F2 = dBCiQ 
and, respectively, for E2 = dB \ Q and Fi = dQ fl B. This contradicts 

(EID-dZSl). 

Therefore, fl must equal B and the theorem is proved. □ 

8. Concluding remarks 

(i) Theorems 12.11 and 13.11 admit a nice probabilistic interpretation, 
e. g., in M^. Informally, it reduces to the following. Consider 
three "Brownian travelers" departing from a point P in the 
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unit ball 5 in R^. The first moves according to the laws of 
standard Brownian motion, the second chooses at random a 
plane through P and follows the Brownian motion in that plane; 
the third chooses at random a line through P and follows the 
Brownian motion on that line. All of these travelers will reach 
the unit sphere S := dB with probability 1. Malmheden's 
theorem asserts that for any portion E G S, the probability 
that the first contact with S occurs in the set E is the same 
for all three travelers, i.e., the observer registering the exiting 
travelers has no way of knowing how they arrived to the unit 
sphere from P. 

(ii) Theorem 12 . 1 1 and . a more general Theorem l3.1l certainly suggest 
connections to integral geometry. Indeed, the Malmheden al- 
gorithm reminds of the inversion formula for Radon transform 
which reconstructs functions from their integrals over hyper- 
planes, or more generally, fc-planes. According to this formula, 
the value of a function at a point coincides, after applying a 
certain power (a half integer in even dimensions) of the Laplace 
operator to the average of the Radon data through the point. 

The Radon inversion formula delivers representation of func- 
tions as continuous sums (i.e., integrals) of so-called plane waves, 
i.e., of functions which are constant on families of parallel planes. 
Similarly, Theorem 13. II can be interpreted in an analogous man- 
ner with the role of plane waves played by " harmonic fc- waves" . 
Here, by harmonic fc- waves we understand functions which are 
harmonic on parallel k- dimensional planes. 

More precisely, fix the natural number k, 1 < k < n. Denote 
by Afc the partial Laplace operator acting only on the first k 
variables: 



For every rotation u G SO{n) we denote by the '"rotated"' 
partial Laplacian: 



It is not hard to show then that the Laplace operator coincides 
with the average of the rotated partial Laplacians, i.e.,: 



k 



dx'j 




n 



k 




SO{n) 
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where doo denotes the normahzed Haar measure on the orthog- 
onal group. Malmheden theorem estabhshes a similar link be- 
tween the solutions of the corresponding Laplace equations: 

Theorem 8.1. Let Q := {x : \x\ < 1} be the unit ball and f G 
C{dQ) and, for every u G SO{n), let u^j denote the (unique) 
solution of the boundary value problem: 



Proof. By construction, the function ou~^, uj G SO{n), is 
harmonic with respect to the first variables xi, ■ ■ ■ , x^, i.e., it is 
harmonic on every cross-section of the ball by a A;-dimensional 
plane parallel to the /c-dimensional plane 11 = {x^+i = • ■ ■ = 
Xn = 0.}. Then, is harmonic on each /c-plane parallel to loW. 
Let X G f2 and let n^j^x be the fc-dimensional plane parallel to 
ujH and passing through x. Since = / on the unit sphere, 
is the harmonic extension of / into ^^^x H ^- When uj runs 
over the whole orthogonal group SO{n), Huj^x runs over all k- 
dimensional planes passing through x and Theorem 3.1 simply 
claims that the average, with respect to all rotations u, of the 
values Uuj{x) equals to the value at x of the harmonic extension 
of the function / into the ball Q, i.e. u{x). This proves fl8.ll) . □ 

The extended Malmheden's theorem for polyharmonic func- 
tions -cf., e.g.. Theorem 16. II for biharmonic functions with k = 1 
- allows a similar interpretation. In this case, the boundary con- 
ditions will be corresponding complete sets of Cauchy data for 
the polyharmonic equation A'^u = 0, i.e., the prescribed values 
of functions and their first — 1 normal derivatives. 

It would be interesting to investigate further whether a simi- 
lar decomposition perhaps holds for other differential operators 
with constant coefficients and rotational symmetry, i.e., some 
operators of the form P{A) where P is a polynomial. 



A^M^(x) = 0, xen, 
Ulu{x) = f{x), X G dQ. 



Then the function 
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(iii) There are various levels at which the converse to the Malmhe- 
den theorem can be formulated. Theorem 14.11 (arguably, the 
most natural one) assumes that Malmheden's procedure pro- 
duces the solution to the Dirichlet problem with arbitrary data 
at all points of a convex domain Q. The conclusion is then 
that n is a ball. Theorem 17. II is an attempt to obtain the same 
conclusion under weaker assumptions: Q is assumed to be star- 
shaped and Malmheden's procedure is only assumed to produce 
the desired results at one point. However, it required an extra 
assumption of central symmetry. There are several other venues 
of interest one may pursue here. To fix the ideas, let f2 be a 
convex domain in M", P E Q, a fixed point. J = Jp, as be- 
fore, is the involution of the boundary F := dQ determined by 
chords through P. As in ^ / = Ip{T) denotes the subspace 
of functions in C(r) invariant under J. We have noted earlier 
in ^ that if P is a Malmheden point of Q (i.e., Malmheden's 
algorithm applied to every / G C{T) produces the value Dpf 
at P of the solution to the Dirichlet problem for the Laplacian 
with data /), then the following hold. 

(a) The two measures on F, dAp = the normalized subtended 
angle from P and Wp := w{-,r,P) := the harmonic mea- 
sure on r at P produce the same results acting on all func- 
tions in /p(r). (Incidentally, this is equivalent to the rela- 
tion between the solid angles of double cones with vertex 
at P and the harmonic measure of the two surface portions 
they cut out of F — cf. §51 Corollary 15.11) . 

(b) For every function / G /p(F) and every linear polynomial 
i, it follows that Dp{e ■ f) = £(P) • Dplf) — cf. Corollary 
15. 2[ Proposition 15.11 

Proposition 15.11 proves the converse, namely: if fl, P are such 
that (a) and (b) hold, then P is an M-point of fl. However, we 
do not know the exact relationship between (a) and (b). 

(iv) Property (a) of a domain Q is on its own somewhat a mys- 
tery and holds, perhaps, the key to a deeper understanding of 
Malmheden's theorem. Theorem 17.11 yields that (a), together 
with the additional hypothesis that Q is symmetric about point 
P, imply that is a ball centered at P. Without this assump- 
tion the conclusion fails. A similar problem in two dimensions 
with the subtended angle measure replaced by a weighted ar- 
clength measure was settled in [8, Thm. 3.29(i)], where it was 
shown that for those problems there are other solutions be- 
sides circles. We must exercise similar caution for our problem 
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as well, since there are star-shaped domains fl (e.g., bicircular 
curves) in the plane for which the subtended angular measure 
from a point A & Q equals to the harmonic measure wb evalu- 
ated at another point B E Q. 
The following example is due to the third author (H. S. Shapiro, 
2006, unpublished). 

Proposition 8.1. There exists a planar domain Q, star-shaped with 
respect to the origin, and a point P 0, P E Q, so that for any subarc 
[/ C r, r := dfl, the angle that U subtends at the origin is equal to 
27rwp{U,Q, P) = 2ttwp{U). (As in previous sections, wp{U) denotes 
the harmonic measure on T evaluated at P.) 

Proof. Consider the polynomial q{z) := az"^ -\-z + a, 0<a< 1/2. 
q is univalent in the unit disk D and maps it conformally onto a 
simply connected domain Q. For z & T, z — e'^, we have q{z) = 
z[a{z + l) + 1] = (1 + 2a cos 6') z. Since 1 2a cos 61 > for all 6, 
q preserves arguments of each z E T. Thus Q is star-shaped and, 
moreover, for any subarc U G T, the angle subtended by U at the 
origin O is equal to the angle that its preimage U' = q^^{U), U G T 
subtends at O. The latter, of course, equals to the length ([/') = 
27rw ([/',D, O) = 2ttwo (U'), the harmonic measure of U' in B evalu- 
ated at the origin. By the conformal invariance of harmonic measure, 
27rw(t/',D,0) = 27Tw{U,n,q{0)) and setting P = q{0) (= a) we are 
done. □ 

B. Gustafsson (a personal communication) has recently character- 
ized all simply connected domains Q C C, where the subtended angle 
measure from an interior point A e Q equals a linear combination of 
the harmonic measures at finitely many other (fixed) interior points. 
But this is the beginning of another tale. 
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